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ABSTRACT 


A  mechanism  Is  proposed  as  a  means  of  passively  leveling  an  optical 
assembly  about  two  horizontal  axes  within  a  tolerance  of  one  arc-second. 
The  mechanism  Is  designed  so  that  It  will  stimulate  a  spherical  pendulum 
for  Input  motions  up  to  100  arc-seconds.  In  order  to  Insure  the  operation 
of  the  mechanism  and  establish  the  proper  geometry,  three  analytical 
models  were  developed.  Planar  and  spatial  models  which  assumed  rigid 
links  with  pin  Joints  and  a  planar  model  assuming  deformable  links  with 
pin  and  clamped  Joints  were  studied.  These  results  are  compared  with 
experimental  data  gathered  for  three  geometry  combinations. 

It  Is  shown  that  the  results  from  the  planar  and  spatial  analytics 
agree  to  within  10%,  with  the  spatial  model  predicting  less  error.  Hence, 
the  less  complicated  planar  analysis  may  be  performed  with  little  loss  In 
accuracy.  Neither  rigid  link  models  compare  well  with  the  experimental 
data,  with  as  much  as  an  80%  discrepancy  In  some  cases.  When 
deformable  links  are  Included  In  the  planar  model,  the  results  are  within 
16%  of  the  experimental  results.  This  Indicates  that  the  resistance  to 
bending  at  the  points  can  not  be  neglected  In  modeling  the  mechanism. 
The  mechanism  to  act  like  a  spherical  pendulum,  and  guidelines  are 


established  for  the  design. 


Fig  1.1a  depicts  the  elements  of  this  concept.  The 
pellet  cerries  e  4  x  4  x  4  matrix  of  detectors  which  ere 
targeted  by  e  laser.  An  entire  row  of  detectors  is  targeted 
simultaneously  by  a  horizontal  thin  plane  of  light.  This 
plene  of  light  is  created  by  sending  the  laser  beem  through 
a  vertically  mounted  cylindrical  lens  of  short  focal  length 
as  In  Figure  1.1b.  An  entire  column  of  detectors  is  tar¬ 
geted  simultaneously  by  rotating  the  eyllndrlcel  lens  into  a 
horizontal  position  as  in  Fig  1.1c. 

Figures  1.2  and  1.3  show  the  emitter  system  in  more 
detail  and  its  passive  leveling  suspension  which  is  Intended 
to  hold  the  Instrument  level  to  within  1  arc  sec. 

1 .2  Variables  in  a  Precision  Two  Axis  Level 

Figure  1.2  is  en  Idealized  suspension  of  flexible  fila¬ 
ments  which  converge  to  a  point.  Such  a  pendulum  will  pro¬ 
duce  a  level  platform  if  the  C.G.  is  directly  below  the  ver¬ 
tex  of  the  system  and  if  all  the  suspension  elements  are 
equal  length  and  have  infinite  flexibility.  In  reality, 
none  of  these  conditions  is  met.  In  the  light  of  the  error 
goal  of  0q  -  1  arc  second,  the  minute  deviations  from  the 
ideal  geometry  can  have  significant  affects.  To  begin  with, 
there  is  no  true  vertex  to  the  suspension  filaments  because 
the  wires  must  be  at  least  .006  in  diameter  to  support  even 
the  low  weight  of  the  lasers  and  lenses.  The  inset  of  Fig¬ 
ure  1.2  demonstrates  that  Instead  of  a  vertex,  the  wires  ere 


actually  separated  by  a  distance  6  so  that  the  suspension  is 
of  the  kinematic  type  shown  by  Fig  1.3.  The  upper  plate 
represents  the  housing  of  the  instrument.  When  the  housing 
Is  tipped  through  an  angle  8^  the  bottom  platform  will  not 
remain  level  but  will  experience  a  tilt  8^  .  Although  6^  is 
smaller  than  8^ ,  It  could  nevertheless  exceed  the  tight 
objective  of  8^  <  1  arc  sec.  The  disruptive  role  of  unequal 
wire  lengths  and  tolerances  in  C.G.  location  will  contribute 
to  the  errors  in  the  Instrument.  The  purpose  of  the  initial 
research  was  to  determine  the  role  of  these  tolerances  to 
the  errors  in  Instrument  suspensions. 

It  is  doubtful  that  the  Instrument  will  be  mounted  on  a 
stationary  surface  since  it  will  be  used  in  a  production 
environment.  It  is  therefore  essential  to  minimize  any 
oscillations  of  the  optics  in  order  to  maintain  the  tight 
constraints  previously  imposed  on  the  system.  So  as  not  to 
interfere  with  the  mechanism,  the  damper  must  be  free  of 
friction.  Figure  1.4  illustrates  an  air  damper  consisting 
of  two  cylinders,  the  inner  attached  to  the  bottom  plate  of 
the  mechanism  while  the  outer  cylinder  is  attached  to 
ground.  Although  simple  in  design,  oscillations  in  both  the 
X  and  Z  directions  can  De  effectively  damped.  Air  damping 
is  more  difficult  to  achieve  due  to  the  low  viscosity  of  the 
air,  while  on  the  other  hand,  it  does  not  require  seals  or 
submersion  of  the  optics  as  would  be  required  with  oil  or 
some  other  medium. 


Spherical  Joints 


Cylindrical 
Lenses  to  Produce 
Two  Perpendicular 
Planes  of  Light 


'2  3\  4 


Cable 


Ternary  Link 


Figure  1.3 


Schematic  representation  of  the  instrument 
indicating  e  3  dimensional  mechanism. 


Cylinder  Attached 
to  Link  5 


Cylinder  Attached 
to  Ground 


Figure  1.4  Cylindrical  air  damper  for  the  two  axes  level 


The  lest  section  of  this  report  focuses  on  the  analysis 
end  experinentetion  required  to  specify  eppropriate  geometr¬ 
ical  relationships  to  facilitate  the  design  of  such  a 


damper 


CHAPTER  2 


PLANAR  ANALYSIS 


2 . 1  Ob  j actives 

The  goal  of  this  analysis  and  those  to  follow  is  to 

determine  the  influence  of  the  geometry  of  the  mechanism  on 

positioning  errors.  The  positioning  error  is  defined  to  be 

a  deviation  of  the  bottom  plate  from  the  horizontal.  This 

angle  will  be  represented  throughout  as  0  .  The  results  of 

o 

this  analysis  will  be  used  to  determine  the  feasibility  of 

this  type  of  mechanism  to  remain  level  and  the  geometric 

constraints  which  must  be  satisfied  in  order  to  minimize  0  . 

o 

2.2  Development  of  Equations 

The  mechanism  under  investigation  is  a  3  dimensional 
linkage.  However,  the  problem  will  be  addressed  initially  by 
reducing  the  general  case  to  2  dimensions.  The  simpler  case 
offers  expedient  insight  into  the  significant  parameters  of 
the  system,  which  for  this  analysis  will  be  the  lengths  of 
the  input  link  (6)  and  output  link  (L),  height  of  the 
mechanism  (h),  the  centering  of  the  optics  (A)  and  tne  dis¬ 
tance  that  the  optics  lie  in  a  vertical  plane  from  Che 


output  link  (e).  The  2  dimensional  mechanism  is  shown  in 
Figure  2.1. 

The  assumptions  used  in  this  analysis  are: 

*  The  cables  are  infiuitly  flexible 

*  Frictionless  pin  joints 

*  Constant  link  lengths 

*  Links  2  &  4  are  two  force  members. 


The  model  is  reduced  to  a  four-bar  linkage  with  tne 
input  rotation  being  link  1,  as  shown  in  Figure  2.2.  A  mass 
is  fixed  to  link  3  at  some  known  location,  G,  relative  to 


p  in 

Joint  A. 

The  kinematics  of  the  linkage  can 

be  written 

as . 

UAC)2  -  (Ax-Cx)2  +  (Ay-Cy)2 

(2.1) 

(1BD)2  -  (Bx-Dx)2  +  (By-Dy)2 

(2.2) 

(Adc)2  "  (Dx-Cx)2  +  (Dy-Cy)2 

(2.3) 

The 

location  of 

the  center  of  gravity  is  fixed  by 

link  5. 

(Acc)2  "  (Gx-Cx)2  +  (Gy-Cy)2 

(2.4) 

(AGd*2  “  (Gx“Dx>2  +  (Gy-Dy)2 

(2.5) 

Referring  to  Figure  2.2,  for  static  equilibrium  the 
forces  and  moments  on  link  3  are  conserved  which  yields: 


E  F  -  0 
3 

fAC%/C  +  '  ■‘T  '  ° 


(2.6) 


s 


Where , 


R  x  t  e  -  R  x  m2  i  m  0 
D/C  BD  B/D  G/C  gJ 


*A/C  "  l 


( Ax- Cx )—  (Ay-Cy )- 

*  1  T  1  J 


AC 


AC 


*B/D 


( Bx- Dx)~  (By- Dy )— 

e  1  »  J 


’BD 


’BD 


-  (Gx-Cx)t  ,  (Gy- Cy )~ 

eG/C  *  1  +  »  J 


GC 


GC 


Insertion  of  the  unit  vectors  into  Equation  2 


f  t 

Ax  -  Cx )  +  -r~— (  Bx-  Dx  )  -  0  :i 
*AC  abd 


Ay-Cy )  +  -j-^(By-Dy)  -  mg  -  0  :j 


AC 


BD 


and  Equation  2.7  becomes: 


-T-—  j(Dx-Cx)(By-Dy)-(Dy-Cy)(Bx-Dx)\  -mg(Gx-Cx)  -  0 
*BD  1  0 


Equations  2.1  -  2.5  and  2.8  -  2.10  describe 
unknowns  for  the  position  solution  of  link  3 
points  A,  B  and  the  geometry,  links  1,  2,  3  and 
solution  was  carried  out  by  solving  the  eight 
equations  for  points  C,D  and  G  numerically  using 


(2.7) 


.6  yields: 

(2.8) 

(2.9) 

(2.10) 

the  eight 
(given  the 
A).  The 
non- 1 inear 
a  Newton- 


Rapson  Method 


The  ouput  eagle  0^  can  Chen  be  found  from  the  tangent  of 
points  D  and  C. 


2.3  Results  from  Planar  Analysis 


In  order  to  generalize  the  results,  the  geometry  will 

be  described  as  ratios  of  link  3  or  "L" .  This  analysis 

showed  that  for  small  angles,  the  input  (0  )  output  (6  ) 

1  o 

relationship  was  linear  as  would  be  expected  tor  small 

e 

o 

motions.  The  ratio  of  has  been  labeled  K  and  is  the 

ei 

measure  of  the  performance  of  the  mechanism. 


The  K  ratio  is  plotted  against  the  ~  ratio  in  Figure 

L 

2.3  for  several  cases  of  j-.  As  can  be  seen,  K  goes  to  zero 
only  when  the  input  link,  (6),  goes  to  zero  (a  true  pendu¬ 
lum)  or  the  height  h  goes  to  infinity  (the  effect  of  6  is 
lessened).  But  there  is  a  physical  limitation  of  the  size 
of  both  6  and  h.  The  construction  of  the  top  link  forces  6 

to  be  at  least  one  wire  diameter.  (See  Insert  on  Figure 

6 

1.2)  Therefore,  the  only  way  to  decrease  —  is  to  increase 

L 

the  length  of  the  output  link.  For  packaging  purposes,  it 

is  desirable  to  minimize  the  height  of  the  mechanism.  From 

Figure  2.3,  it  is  therefore  advisable  to  generate  a  mechan- 

h  h 

ism  which  has  an  —  ratio  about  unity  since  —  much  greater 

L  la 

than  unity  offers  little  Improvement  on  K. 


The  effect  of  is  illustrated  in  Figure  2.4,  where  it 
A 

can  be  seen  that  —  may  vary  up  50%  without  introducing  posl- 


0.002  0.004  0.006  0.008  0.010 


Figure  2. 


6 

c 


K  vs  6/L  as  a  function  of  h/L 
from  the  planar  analysis. 


tioning  errors.  At  first  this  result  seemed  spurious,  but  i 

look  et  the  conf igure tioo  of  the  mechanism  with  6  «*  0 

o 

explains  the  results.  The  inset  of  Figure  2.4  shows  that 
during  assembly  the  instant  center  of  link  1  with  respect  to 
link  3  is  brought  to  a  point  directly  above  tne  mass  center 
of  the  mechanism,  while  the  wire  lengths  are  altered  such 
that  the  bottom  link  is  horizontal.  By  a  simple  tri¬ 
gonometric  analysis,  it  can  be  shown  that  this  assembly  pro¬ 
cedure  causes  the  distance  from  the  center  of  mass  to  the 
instant  center  to  be  independent  of  A.  (The  distance 
between  links  1  and  3  is  fixed  by  h.) 

Since  the  angular  velocity  ratio  tor  any  two  links  in  a 
mechanism  is  Inversely  proportional  to  the  distances  from 
the  Instant  centers  in  the  frame  about  which  the  links  are 
rotating  to  the  Instant  center  which  is  common  to  the  two 
links,  the  assembly  procedure  fixes  the  velocity  ratio.  The 
angular  velocity  ratio  is  simply  K  since  the  relative  rota¬ 
tions  are  small  in  this  mechanism.  Thus,  K  remains  unal¬ 
tered  because  the  distances  from  the  instant  center  to  links 
1  and  3  are  independent  of  A  and  therefore  fixing  the  velo¬ 
city  ratio,  K,  by  this  constant  distance  ratio. 

e 

The  effect  of  ~  on  K  is  shown  in  Figure  2.5.  It  is 
Important  to  note  that  the  mass  center  should  lie  below  the 
bottom  link,  since  K  increases  as  the  output  link  length  e 


decreases 


With  a  negative  value  of  e,  K  Increases  monoton- 


ically,  whereas  a  large  positive  value  for  e  significantly 
reduces  K.  Since  the  value  of  —  is  a  design  choice,  it  is 
recommended  to  make  this  value  as  large  as  possible. 


T!*P 
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CHAPTER  3 


SPATIAL  ANALYSIS 


3.1  Ob  j  ectlvea 

While  the  2  dimensional  analysis  is  a  good  indicator  of 
the  influence  of  the  geometry  on  the  output  ,  a  3  dimen¬ 
sional  analysis  will  show  the  rotations  out  of  the  plane  of 
the  2  dimensional  mechanism.  In  this  section  we  will  focus 
on  the  static  equilibrium  position  of  the  mechanism  in  3 
dimensions  subject  to  input  rotations  (0^)  about  the  X  and  Z 
axis  as  Illustrated  in  Figure  3.1.  The  results  will  then  be 
compared  to  the  results  from  the  previous  section  in  order 
to  understand  the  validity  of  the  two-dimensional  analysis. 

3.2  Development  of  Equations 

Multiple  starts  were  made  to  this  problem  in  order  to 
find  an  approach  with  sufficient  accuracy  to  indicate  rota¬ 
tions  on  the  order  of  arc  seconds,  while  simultaneously 
holding  the  computational  time  to  a  minimum.  The  first 
attempt  was  to  use  IMP  [1],  a  general  purpose  mechanism 
analysis  computer  program,  but  the  solution  was  too  coarse 
Tl ]  She th ,  P.N.  and  Vicker,  J.J.  Jr.,  "IMP  (Integrated 


,•  VV  .  •  ,  *  .s  •-  v  s’  .  v. 


to  b«  of  any  value.  The  vector  loop  approach  producer  cou¬ 
pled,  transcendental  equations.  Uhlle  this  creates  no  diffl- 
cultles  In  planar  analysis,  standard  numerical  techniques 
experience  difficulties  in  converging  on  a  solution  to  the 
large  number  of  equations  required  for  spatial  analysis. 
The  situation  does  not  change  when  small  rotations  are 
assumed  about  the  equilibrium  position  because  the  equations 
remain  linearly  coupled  by  cross  products.  Since  the  vari¬ 
ables  are  small,  again  convergence  problems  arise  In  numeri¬ 
cal  equation  solvers. 


Minimization  of  potential  energy  was  explored  as  an 
alternative  means  of  developing  the  equations  of  motion. 
The  creation  of  multiple  partial  differentials  which  could 
not  be  determined  algebraically  forced  the  use  of  Lagrange 
Multipliers  to  evaluate  the  partlals.  When  the  Lagrange 
Multipliers  are  coupled  with  the  constraint  equations,  the 
number  of  equations  Is  Increased  to  30. 


The  final  approach  contained  herein  Is  less  elegant  In 
nature  but  requires  significantly  fewer  constraint  equations 
(15).  These  equations  are  still  non-linear,  but  the  equa¬ 
tion  solver  converges  rapidly  because  the  variables,  the 
coordinates  of  the  spherical  joints,  are  relatively  large. 


Mechanisms  Program),  A  Computer-Aided  Design  Analysis 
System  for  Mechanisms  and  Linkages",  Journal  Engineering 
for  Industry,  Trans.  ASHE,  Ser.  B,  Vol.  99,  no.  3, 
August  1977,  pp.  773-778. 
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This  approach  to  the  spatial  analysis  is  basically  the 
sane  as  the  planar  analysis,  except  It  allows  motion  In 
three  directions;  that  Is,  the  mass  ot  the  optics  Is  at  some 
point  G  on  the  ternary  output  link  (as  shown  In  Figure  3.1) 
and  the  link  lengths  are  constant  with  frictionless  spheri¬ 
cal  joints.  The  input  rotation  Is  about  the  X  and  Z  axis 
located  at  the  center  ot  the  input  link.  Rotation  ot  the 
input  link  about  the  Y  axis  does  not  alter  the  magnitude  ot 
the  output,  In  that  it  merely  rotates  the  coordinate  system. 
The  analysis  is  predicated  on  knowing  the  initial  coordi¬ 
nates  of  points  A,  B  and  C.  Following  any  input  (9^),  the 
resulting  coordinates  of  A,B  and  C  are  calculated  via  a 
rotation  matrix. 

The  first  equations  ot  constraint  are  those  ot  constant 


link  lengths.  Figure  3.1. 

2  2  2  2 
(1  )  *  (Dx-Fx)  +(Dy-Fy)  +(Dz-Fz)  (3.1) 

D  r 

UDE)2  -  (Dx-Ex)2+(Dy-Ey)2+(Dz-Ez)2  (3.2) 

UEp)2  -  (Ex-Fx)2+(Ey-Fy)2+(Ez-Fz)2  (3.3) 

(*AD)2  -  (Ax-Dx) 2+(Ay- Dy) 2+(Az-Dz) 2  (3.4) 

UBE)2  -  (Bx-Ex)2+(By-Ey)2+(Bz-Ez)2  (3.5) 

UCF)2  "  <Cx-Fx)2+(Cy-Fy)2+(Cz-Fz)2  (3.6) 

(AGD)2  "  (Gx-Dx) 2+(Gy-Dy) 2+(Gz-Dz) 2  (3.7) 

UGE)2  -  (Gx-Ex)2  +  (Gy-Ey)2  +  (Gz-Ez)2  (3.8) 

UGF>2  -  (Gx-Fx)2+(Gy-Fy)2+(Gz-Fz)2 


(3.9) 


For  static  equilibrium,  Figure  3.2 


E  F$  -  0 

fAD*A/D  +  fBE*B/E  +  fCF*C/F  " 


(3.10) 


^D/G  *  fAD*A/D  +  ^E/G  X  fBE*B/E 
+  \/G  X  fCF*C/F  "  ° 


(3.11) 


When  the  unit  vectors  are  introduced  into  Equation  3.10  and 
the  multiplications  carried  out,  these  equations  become: 


£  f  f 

7^(Ax-Dx)+7H(Bx-Ex)+7H(Cx-Fx)  -0  :  i 

XAD  XBE  XCF 

y^(  Ay- Dy)  l--B-(  By- Ey  )+■?--( Cy-Fy)  -mg  :J 
*AD  XBE  *CF 

f  £  f  __ 

*j“—(  Az-D*)+-— ( B  z- Ez)  +«j— —  ( Cz- F  z)  -0  :k 
AD  *BE  XCF 


(3.12) 


(3.13) 


(3.14) 


Taking  the  cross  products  in  Equation  3.11  yields  the  foi- 
1  owing : 

£ 

T~1  (Dy-Gy)(Az-Dz)  -  (Dz-Gz)  (Ay-Dy)l 
*AD  } 

£ 

+  •j——  $(  Ey-  Gy )  (Bz-Ez)  -  (  Ez- Gz  )  (  By- Ey 
BE  1  ’ 

£  _ 

+  ^(Fy-Gy)  (Cz-Fz)  -  ( F z- Gz )  ( Cy- Fy )|  -0  :i  (3.15) 


— ^  J(Dz-Gz)  (Ax-Dx)  -  (Dx-Gx)  (Az-Dz)? 

*  AD  ^ 

+  T~“  l(Ez-Gz)  (  Bx-  Ex  )  -  (  Ex-Gx  )  (  B  z- Ez  )\ 

*BE  3 

+  T”^?I’GzHCx-Fx)  -  (Fx-Gx)  (Cz-Fz)^  -0 


(3.16) 


-^£{(Dx-Gx)  (Ay-Dy)  -  (Dy-Gy)  (Ax-Dx)^ 

XAD 

+  T^^Ex-Gx>^By-Gy)  “  (Ey-Gy ) (Bx-Exft 
XBE 

+  j^-^(Fx-Gx)  (Cy-Fy)  -  (Fy-Gy) ( Cx-Fx)^ -0  :k 


(3.17) 


Equations  3. 1-3.9  and  3.12-3.17  are  then  solved  numeri¬ 
cally  for  the  coordinates  of  points  D,  E  and  F  on  the  output 
link.  The  output  rotations  are  then  calculated  using  Rodri¬ 
gues"  formula  [2]  to  find  the  rotation  of  the  normal  to  the 
output  link  from  the  Initial  assembly  position.  Rodrigues" 
formula  is  used  rather  than  a  general  displacement  matrix  to 
reduce  computational  time. 


3.3  Comparison  of  the  Spatial  and  Planar  Results 


First,  Rodrigues"  formula  is  only  valid  when  a  body  is 
rotating  about  a  fixed  axis  without  any  translations.  This 
was  verified  oy  calculating  the  displacements  of  the  output 
link  using  a  general  displacement  matrix,  and  then  noting 
that  the  displacement  terms  were  equal  to  zero. 


The  primary  results  of  this  analysis  are  displayed  in 
Figure  3.3,  where  K"  is  plotted  against  ~  for  different 
values  of  The  ratio  K"  is  the  magnitude  of  the  rotation 

L 

of  the  output  link  about  the  X  and  Z  axes  over  the  magnitude 
of  the  input  rotation.  The  3  dimensional  analysis  of  the 

l 2]  Bisshopp ,  K.E.,  "Finite  Rotations  of  a  Rigid  Body  and 
the  Airplane  Landing  Gear  Problem",  Jnl.  Mechanisms, 
Vol.  3,  Pergamon  Press,  1968  ,  pp.  203-208  . 


c 


Figure  3.3 


K'  vs  6/L  as  a  tunction  of  h/L  from 
the  spatial  analysis  with  comparisons 
to  the  planar  analysis. 


mechanism  Indicates  a  smaller  K  than  that  predicted  by  the  2 
dimensional  analysis  due  to  rotations  on  the  order  of  hun¬ 
dreds  of  arc-seconds  of  the  output  link  about  the  Y  axis. 

Figure  3.4  illustrates  K"  verses  for  various  £ 

ratios,  and  a  comparison  with  Figure  2.5.  As  previously 

shown  in  Figure  3.3,  the  spatial  analysis  yields  a  lower  K 

than  does  the  planar.  It  can  be  seen  that  the  difference 

between  the  spatial  and  planar,  at  each  is  similar  to 

that  seen  in  Figure  3.3,  and  that  the  curves  tor  the  spatial 

h  z 

and  planar  models  begin  to  coalesce  as  ~  or  —  go  to  infin¬ 
ity. 

6 

Since  the  two  models  predict  the  same  performance  as  ~ 

la 

h 

goes  to  zero  or  —  goes  to  infinity,  we  can  conclude  that  the 

la 

planar  analysis  is  a  suitable  substitute  for  the  spatial 
model,  with  the  understanding  that  K  found  from  the  planar 
is  larger  than  K'  found  from  the  spatial  analysis. 

Figure  3.5  illustrates  the  change  in  one  of  the  output 
link  lengths,  while  the  two  others  are  held  equal.  The  fig¬ 
ure  shows  that  K'  actually  decreases  below  the  nominal 
value,  all  lengths  (L)  equal,  as  one  length  is  increased  to 
roughly  1%  larger  than  the  others.  Let  the  input  rotation 
be  about  the  axis,  illustrated  in  the  inset  of  Figure  3.5, 
which  is  measured  from  the  positive  X  axis  parallel  to  the 
link  which  is  changing,  described  by  the  angle  a.  The  lower 
curve,  for  which  the  change  in  K'  is  very  severe,  is  the 


010 


analysis 


Input  rotation  at  a 


90  degrees  from  the  length  that  das 


been  changed.  When  the  rotation  is  parallel  to  the  length, 
o*0  degrees,  the  minimum  is  the  nominal  value. 


By  looking  down  the  a  *  0  degree  axis  in  Figure  3.5, 
one  can  see  that  the  apparent  link  lengths  are  not  altered 
by  the  change  in  length,  end  the  mechanism  performs  as 
though  no  change  had  been  made.  On  the  other  band,  looking 


down  the  90  degree  axis,  one  sees  the  full  effect  of  the 

change  in  length,  and  hence  a  decrease  in  K".  The  minimum 

AL 

value  of  K" ,  at  ~ —  *  1%,  is  less  than  K'  if  all  links  are 
increased  by  1%.  This  perhaps  could  be  advantageous  if  the 


mechanism  was  used  to  level  the  optics  about  only  one  axis. 


The  results  of  the  spatial  analysis  are  very  similar  to 
those  of  the  planar  analysis.  The  planar  analysis  is  much 
shorter,  and  the  results  equally  as  informative  as  the  spa¬ 


tial  analysis.  The  effects  of  one  link  length  changing 


shows  that  K'  decreases  below  the  nominal  minimum 


CHAPTER  4 


PLANAR  MODEL  ACCOUNTING  FOR  WIRE  STIFFNESS 


4 . 1  Obj  ec tivea 

As  previously  mentioned,  links  2  &  4  are  made  of  cable 
with  less  than  ideal  joints  at  the  ends.  They  have  been 
assumed  to  have  infinite  flexibility  so  that  the  bending  in 
the  cables  occurs  at  an  infinitesimal  area  at  the  the  con¬ 
nection  to  the  plates. 


For  this  section  the  cables  trill  have  some  finite 
stiffness  and  are  subjected  to  shearing  forces  and  moments 
applied  to  the  ends.  A  planar  model  of  the  mechanism  is 
developed  here  which  includes  the  stiffness  in  the  cables  as 
a  system  parameter.  It  will  be  shown  that  Including  this 
effect  in  the  kinematic  representation  of  the  mechanism 
leads  to  results  which  will  agree  well  with  the  experimental 


data  gathered  in  the  next  section 


Development  of  Equations 


The  following  easumptioi 


*  Negligible  strain  energy  Is  stored  in  the  cables  at  the 
assembly  position. 


*  The  cables  are  Euler-Bemoulli  beams  such  that  the  cur¬ 
vature  is  proportional  to  the  bending  moment. 


*  Deflection  of  the  beams  Is  small  such  that  bending  does 
not  alter  the  length  of  the  cables. 


*  Strains  In  the  Individual  wires  of  the  cable  are  small 
and  the  axial  and  bending  responses  are  additive. 


Friction  and  twisting  in  the  cables  are  neglected, 


4.2.1  Force  and  position  solution 


Upon  assembly  of  the  mechanism,  the  cables  are  plasti¬ 
cally  deformed  around  the  junction  at  the  input  link  and 
subsequently  assume  an  unstressed  kinked  shape  at  this  point 
if  disassembled.  The  assumption  of  no  Initial  stored  strain 
energy  Is  then  an  approximation  to  the  state  of  the  cables 
at  assembly.  This  assumption  will  lead  to  a  pinned  boundary 
condition  at  the  input  link  whereas  a  clamped  boundary  con¬ 
dition  will  be  used  at  the  output  link.  The  cables  will  be 
constrained  to  make  small  deflections  about  the  assembly 


position  In  order  to  simplify  the  analysis.  These  assump¬ 
tions  lead  to  a  schematic  representation  of  the  mechanism  as 
illustrated  in  Figure  4.1. 


Let  G.  end  9  be  measured  from  the  Initial  configura- 
A  B 


tion  of  the  mechanism  upon  assembly.  0  and  6  must  be 

A  B 


altered  by  the  rotation  of  the  Input  link,  6  ,  In  order  to 


maintain  the  axial  force  parallel  and  the  shear  force  per¬ 
pendicular  to  the  undeformed  cable.  The  subsequent  angles 
after  a  rotation  of  9^  will  be  4^  and  4g  as  shown  In  Figure 


4.2. 


4  -  0  -  0 

A  A  1 


4-9+0 
B  B  1 


For  static  equilibrium,  the  force  and  moment  balance  on 
the  output  link  yields. 


E  F  -  0 
3 


HACO**A  +  PASin*A  +  V°S*B  "  PB*la*B  "  °  <4*L) 


PACO-*a  “  HA*l“*A  +  PBC°*^B  +  HB*ln*B  "  "  0  (4*2> 


E  «a  -  ° 

•"a  •  \  +  VC<,,<W  +  “»L“"(  W 


L 

-mg  (T  -  6)cos0  +  eslnG  -  0  (4.3) 

Z  o  o 

Referring  to  Figure  4.3,  the  deflections  of  the  cables 
are  along  a  line  perpendicular  to  the  undeformed  initial 
position,  so  that  point  A  can  be  calculated  from 


A 


A. slnd . 


W.COS  + 


Similarly, 


\  ■  -V°**A  +  "A*1"** 


B  “  6cos9  +  Jl.slnA  -  w  coil 
x  1  B  B  B  B 


By  -  6smei  -  Xbco»*b  -  wBam^B 


The  length  of  the  notton  plate  remains  unchanged,  thereby 


connecting  the  deflections  w  and  w  . 

A  D 


2  2  2 
(B  -A  )  +  (B  -A  )  -  X 

xx  yy  o 


(4.4) 


B  -  A 

— ^ - -X-  »  tanQ 

B  -  A  o 


(4.5) 


4.2.2  Derivation  of  beam  equa tl 


ons 


The  cables  have  been  assumed  to  be  beams  which  are 
loaded  primarily  In  tension.  The  Important  equations  will 
be  derived  for  one  cable  with  the  second  cable  being  loaded 
In  a  similar  way. 


Referring  to  Figure  4.4,  and  talcing  moments  about  point 


P  * 


M  -  M  +  P(w  -  w)  -  H  (A  -  u)  *•  0 


Y 


Figure  4.3  Deflection  diagram  for  both  cables  as  used 

for  the  position  solution  of  points  A  and  B. 


Sine*  the  cables  hive  been  assumed  to  be  Euler  Bernoulli 


beams  with  small  curvature,  we  can  then  write: 


2 

d  w 


M 

El 


d  M„  2  2  P 

- - —  -  p  M  ■  0  where  P  »  —  (4.6) 

,  2  u  El 

du 

where  El  Is  the  flexural  rigidity  of  the  cable.  The  general 
solution  of  Equation  4.6  is 

M  ■  C^coshpu  +  C^slnhpu 
o  r 

d2w  82 

- j  »  --|~(CLcoshPu  +  C^slnhpu)  (4.7) 

d  u 

As  previously  discussed  in  section  3.2.1,  the  boundary 
conditions  will  make  use  of  the  fact  that  the  cables  are 
nearly  pinned  at  the  input  link  and  nearly  clamped  at  the 
output  link.  The  boundary  conditions  can  be  visualized  with 
the  aid  of  Figure  4.5  where  the  cable  and  plates  are  magnl> 
fled  to  show  the  connections.  The  cable  is  firmly  clamped 
at  the  bottom  plate  via  a  screw  and  washer.  The  cable  at 
the  top  is  drawn  around  a  knife  edge  yielding  a  joint  which 
baa  the  aspect  of  being  sensitive  to  plate  rotation  while 
not  exerting  a  moment  on  the  plate.  So  a  pin  joint  is  a 
reasonable  approximation  to  the  actual  boundary  conditions 


The  cables  are  not  securely  fastened  to 
the  input  link  at  the  knife  edge. 

Since  the  attachment  is  away  from  the 
knife  edge,  a  pinned  condition  occurs. 


The  cables  are  securely  clamped  to 
the  output  link  by  screws. 


Figure  4.5 


Sketch  of  the  instrument  Illustrating 
the  cable  boundary  conditions. 


at  the  top  link.  The  boundary  conditions  can  then  be  stated 


u  -  0 


u  -  Jl 


w  ■  0 


12  -  8 
d  u  o 


The  above  boundary  conditions  Lead  to  the  evaluation  of 
the  constants  In  Equation  4.7. 


C  -  c  -  0 
1  4 


C  -  ~ _ 

2  slnhf3A 


r  ■  a  4-  PMcoshP* 

3  ao  *  PslnhpA 

Then  the  displacement  and  shear  at  the  bottom  of  the  cable 
can  be  found  from  equation  4.7  and  its  derivatives. 


(C^coshp*  +  s lnhpi)  + 


Inserting  the  constants  yields 


From  the  equation  for  the  shear  on  the  cable 


H  d  w 

El  3 

d  u 


(4.8) 


•  %  -  *  .**  , 


v'v*  »*'  .** 


■.'V.v'Va 


11 

p 


1* 

El 


•  -f<c  ^slntipjl  + 


C2  coshf3Jl) 


or 


i 

—t—  -  C.slnnpjl  +  C,  cosh{3A 
P  1  1 


Inserting  the  constants. 


cosnBl 

H  -  -M.p - ~ 

A  JTslnhpA 


(4.9) 


Equations  4.8  and  4.9  are  then  written  for  each  cable. 
With  the  inclusion  of  equations  4.1  -  4.5  there  are  enough 
constraint  equations  to  solve  for  the  equilibrium  position 
of  the  output  link  numerically. 


4 . 3  Results  from  stiffness  analysis 


In  order  to  evaluate  the  magnitude  of  the  stiffness 

effects  on  K,  the  moment  of  inertia  (I)  and  the  modulus  of 

elasticity  (E)  were  found  for  a  0.006  in.  diameter  7X1 

stainless  steel  cable.  The  moment  of  inertia  for  the  cable 

-  1 1  4 

was  calculated  to  be  4.494  *  10  in  ,  while  the  modulus  of 

6 

elasticity  was  determined  to  be  approximately  14  *  10  psl 

6 

[ 3 ]  [ 4 ] •  Stainless  steel  has  a  modulus  of  28-29  x  10  psl, 

but  a  correction  has  been  made  for  slippage  and  the  reduced 

area  of  the  cable  as  opposed  to  a  solid  wire.  These  values 

are  only  approximate  as  no  experimental  results  were 

TTj  Cos  te 1 lo ,  G.A. ,etal," Effective  Modulus  of  Twisted  Wire 
Cables",  Journal  of  the  Engineering  Mechanics  Division, 
ASCE,  Vo  1 .  102,  Ho.  EMI,  February,  1976  ,  pp.  171-181. 

(4]  Shlgley,  J.E.  ," Mechanical  Engineering  Design",  Third 
Edition.  McGraw-Hill,  New  York,  1977. 


obtained  to  confirm  the  moment  of  inertia  or  modulus  of 


elasticity  of  the  cable. 


Figure  4.6  illustrates  the  Influence  of  the  stiffness 

1  6  h  h 

on  ~  as  a  function  of  ~r  and  — .  The  effect  of  —  is  shown  by 

lx  L  L  Li 

successive  plots  starting  with  -  1/2  (Figure  4.6a)  and 

L 

increasing  up  to  j-  *  10,  (Figure  4.6d).  As  with  the  earlier 

la 

h 

analysis,  ~  >  1  has  little  effect  on  the  input  output  ratio 
K.  In  fact,  there  is  no  difference  between  Figures  4.6c  & 


The  ~  ratio  is  plotted  on  each  figure  tor  what  kinemat- 
Li 

6 

ically  would  be  a  pure  pendulum,  ■  0 .  K.  does  not  go  to 

zero  as  it  would  in  the  kinematic  analyses,  but  instead  is  a 

strong  function  of  El.  It  is  apparent  from  Figure  4.6  that 

6 

the  kinematic  influences  at  small  —  ratios  are  negligible 

6 

compared  to  the  stiffness  influences.  As  ~  Increases  the 
kinematic  errors  begin  to  out  weigh  the  stiffness  errors. 

Therefore  the  earlier  analyses  are  valid  for  small  values  of 

6 

El  or  for  large  values  of  -jj. 

In  Figure  4.7  the  effects  of  j-  are  illustrated.  Fig- 

la 

-  „  h  6 

urea  4.7a  through  4*7d  are  Dased  on  “  ■  1  while  y  is  lncre- 

L  la 

mented  on  each  plot.  Again  ^  has  a  strong  Influence  on  K 

la 

and  and  to  minimize  the  positioning  error  of  the  mechanism  a 


large  —  is  recommended. 
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Figure  4.6  1/K  vs  mg*2/EI  it  e 

function 

of  6/L 

for  varying  h/L  from 

planar  analysis 

including  stiffness. 

In  Chapter  2  it  was  shown  kinematically  that  the  CG 
could  be  located  anywhere  up  to  50%  of  L  along  the  output 


link  without  an  associated  error.  The  inclusion  of  stiff* 
ness  into  the  kinematic  model  shows  that  A  may  vary  widely 
without  introducing  any  errors  and  also  that  the  maximum 
limits  are  dependent  on  the  cable  stiffness. 


A 

Figures  4.8a, b  illustrates  ^  of  30%  and  20%  respec¬ 
tively.  The  curves  do  not  differ  from  that  of  ~  -  0,  lndi- 

L 

eating  that  miscenterlng  of  the  CG  is  not  as  critical  as 

indicated  in  Chapter  2.  Because  the  bending  of  the  cables 

is  creating  a  reactive  moment  on  the  output  link  the  sign  of 

A  is  important.  If  A  is  increased  in  the  direction  of  the 

deflection,  (l.e.  to  the  left  in  Figure  4.4)  the  maximum 
A 

value  for  —  is  approximately  30%.  On  the  other  band,  if  A  is 

negative  with  respect  to  the  displacements,  the  maximum  is 

close  to  40%  before  the  cables  are  loaded  in  compression. 

These  results  are  dependent  on  El  and  are  stated  for 
2 


1000,  or  the  first  point  plotted  in  Figure  4.8.  As 


El  Increases  these  values  would  Increase  to  50% 


CHAPTER  5 


EXPERIMENT 


5 . 1  Objectives 

Experimental  results  were  needed  to  Identify  problems 
with  the  mechanism  and  to  give  a  bench  mark  tor  the  analyt¬ 
ics.  Three  mechanisms  were  tested,  two  with  a  large  Input 
range  and  one  for  much  smaller  Inputs.  It  was  desired  that 
these  testa  indicate  the  roles  of  the  geometry  of  the 
mechanism,  hence  three  different  geometries. 

5 . 2  Description  of  Experiments 

The  experimental  procedure  was  divided  into  a  small 
input  test  (for  input  motions  from  10  to  300  arc-seconds) 
and  a  large  input  test  (for  rotations  over  300  arc-seconds). 

The  small  Input  test  is  shown  in  Figure  5.1a.  An  auto¬ 
collimator  was  used  to  monitor  both  input  (0^)  *ad  output 

(0  )  by  use  of  a  beam  splitter.  The  autocollimator  enabled 
o 

a  0.5  arc-second  accuracy  to  be  held  on  both  the  input  rota¬ 


tions  and  the  output  rotations 


Since  Che  autocollimator  bad  only  a  ±  5  arc-minute 


range,  the  input  rotations  during  the  coarse  experiment  were 
monitored  by  using  a  spirit  level  witn  5  arc-minute  gradua¬ 
tions  as  shown  in  Figure  5.1b  Experimental  procedures  con¬ 
sisted  of  rotating  the  lab  stand  to  a  predetermined  0^  via  a 
micrometer  and  then  reading  the  output  angles  by  use  of  the 
autocollimator. 

5.3  Results  from  Static  Positioning  Experiments 

Tests  were  conducted  on  three  mechanisms  which  are 
described  by  the  table  (Table  5.1)  of  non-d lmens ionles s  link 
lengths.  The  major  differences  between  the  mechanisms  were 
the  parameters  (6/L).  Both  of  the  cases  A  and  B  have  rela¬ 
tively  small  input  links  (6)  compared  to  case  C.  The  height 
(h)  of  A  and  B  was  again  much  smaller  than  C  as  indicated  by 
the  h/L  ratio  in  the  second  column.  Plus,  the  stiffness 
mgi.2 

parameter  -g  ■ —  for  case  A  was  much  smaller  than-  cases  B  and 
E  X 

C  due  primarily  to  an  Increased  mass  of  these  mechanisms. 

Table  5.1 

Non-Dimensional  Ratios 
Describing  Experimental  Mechanisms 


Case 

6 

L 

h 

L 

A 

L 

e 

L 

A 

0.0034 

0.722 

0.0 

0.144 

9342 

B 

0.0021 

0.650 

0.0 

0.072 

12160 

C 

0.0197 

1.559 

0.0 

0 . 144 

12160 

5.3.1  Sm«Ll  input  test 


Eight  experioentsl  states  ere  plotted  in  Figure  5.2  for 
the  Case  A  mechanism  along  with  the  theoretical  values  from 
the  planar  kinematic  and  planar  stiffness  models.  At  each 
input  state,  an  error  band  Indicates  the  range  of  the  meas¬ 
urements  while  the  circle  indicates  the  average.  Linear 
regression  yields  the  best  fit  line  through  the  data  and  the 

slope  of  that  line  is  K  .  Comparisons  of  K  to  the 

exp  r  exp 

kinematic  model  of  Chapter  1.  K  and  the  kinematic  model 

kin 

including  cable  stiffnesses,  K  are  then  shown  in  Table 

atif 


5.2. 


It  is  apparent  from  Table  5.2  that  the  kinematic 

analysis  yields  a  K  which  is  extremely  low  in  comparison  to 

the  experimental,  with  a  difference  of  78. 6Z  from  K,  _  to 

k  la 

K  .  The  second  analysis  gives  a  much  better  Indicator  of 
exp 

the  performance  of  the  mechanism  with  only  a  13.4%  differ¬ 
ence  between  K  and  K  , The  primary  reason  for  the 

exp  stif  *  7 

large  change  between  K,  ,  and  K  is  due  to  the  small 

kin  stif 

value  of  .  As  seen  in  Chapter  4  Figure  4.6,  the  effects  of 

L 

6 

the  stiffness  become  predominant  as  —  becomes  small. 


test  for 


Table  5.2 

Comparison  of  Results 
for  Experimental  Mechanisms 


Case 

K 

K, 

%  Diff 

K 

%  Diff 

exp 

kin 

s  tif 

A 

0.0187 

0.0040 

78.6% 

0.0162 

13.4% 

B 

0.0124 

0.0020 

85.0% 

0.0144 

-16.1% 

C 

0.0194 

0.0185 

5.6% 

0.0221 

-13.9% 

(%  Diff.  defined  as  — L-£  x  100%  ) 

Kexp 


5.3.2  Large  input  test 

The  results  from  the  large  input  test  for  case  B  are 

illustrated  by  the  lower  set  of  curves  in  Figure  5.3.  Case 

B  consistently  has  a  smaller  angular  error  ( 0Q )  than  does 

Case  C  as  predicted  by  the  planar  kinematic  analysis  of 

each.  The  kinematic  analysis  again  predicts  a  lower  K  than 

that  found  experimentally,  as  indicated  in  Table  5.2  with  a 

85.6%  difference  of  K,  to  K  .  Although  the  stiffness 

kin  exp 

model  reflects  the  experimental  data  better,  with  a  16.1% 

difference  from  K  ,  to  K  ,  it  also  predicts  a  larger 

stiff  exp 

value  for  K  than  that  measured.  This  is  surprising  in  view 

of  the  case  A  results  where  K,  ,  was  larger  than  K 

kin  exp 

Turning  to  case  C,  illustrated  by  the  upper  set  of 
curves  in  Figure  5.3,  it  is  apparent  that  the  planar 


kinematic  analysis  for  this  mechanism  accurately  predicts 
the  experimental  results  with  a  difference  of  only  5.6%  as 


es  t  for 


shown  In  Table  5.2.  This  shows  that  the  kinematic  errors 

6 

due  to  the  large  ~  are  outweighing  the  errors  due  to  stiff- 

la 

ness.  The  inclusion  of  stiffness  in  the  analysis  increases 

K  to  0.0221,  for  a  difference  of  13.9%  between  K,  and 

stiff  kin 

K  with  the  stiffness  model  again  yielding  the  larger  K. 

s  1 1  f 

5.3.3  Discuss  ion  o  f  results  and  error  sources 

As  mentioned  previously,  the  symbols  on  the  experimen¬ 
tal  curves  reflect  the  average  of  the  data  at  that  point, 
while  the  error  bands  extend  above  and  below  these.  The 
error  band  was  as  much  as  ±  30%,  15%,  and  10%  of  the  average 
values  for  cases  A,  B  and  C  respectively.  For  case  A  the 
large  spread  is  partially  due  to  the  resolution  of  the  auto¬ 
collimator  (±  .5  arc  seconds)  while  a  1  arc  second  measure¬ 
ment  was  being  made.  Hence,  an  accuracy  problem  exists. 
For  cases  B  and  C  the  spirit  level  used  to  monitor  the  input 
could  be  read  to  ±  1  arc  minutes,  while  the  magnitude  of  the 
input  angles  were  on  the  order  of  1  degree.  So  it  can  be 
concluded  that  cases  B  and  C  were  run  with  reasonable  accu¬ 
racy,  while  the  data  for  case  A  could  be  different  from  the 
nominal  values  plotted  in  Figure  5.2  or  listed  in  Table  5.2. 

Since  the  stiffness  analysis  predicts  K  values  greater 
than  the  experimental  for  cases  B  and  C,  it  indicates  that 
the  analytics  are  Incorrect  or  that  some  variable  or  values 
were  Incorrect.  As  mentioned  m  Section  3.3,  the  value  of 
El  was  merely  an  estimate.  If  El  were  decreased  by  a  small 


14  X  10  is 


amount  K  would  decrease  sizably  since 


seen  to 

lie  on  the 

rapidly 

increasing 

portion 

of 

the 

curve 

s  hown 

in  Figure 

3.6b. 

From  this 

figure 

one 

sees 

that  a 

decrease  of  E  from 

14  to  12 

6 

X  10  ,  a 

decrease 

of 

14%, 

would 

change 

K  approximately  6%. 

Also,  the 

difference 

between  the 

planar  and  spatial  analyses  from  Chapter  3  was  roughly  10%, 
thereby  reconciling  the  differences  between  the  experimental 
results  and  the  results  from  the  kinematic  model  Includ lng 
stiffness  as  being  dependent  on  the  value  tor  E  and  the 
rotations  out  of  the  plane  of  the  planar  analysis. 

It  should  be  noted  that  errors  could  be  introduced  in 
the  planar  analysis  Including  stiffness  where  small  deflec¬ 
tions  were  assumed.  Inherent  In  this  assumption  is  that  the 
distance  along  the  undeformed  cable,  (u)  from  Chapter  4  Fig¬ 
ure  4.7,  remains  constant  and  the  deflections  (w)  are  from 
the  undeformed  cable.  The  actual  length  of  the  cable  Is 
2  2 

X  •  \|u  +  w  .  While  this  is  a  small  effect,  it  would 

change  K  by  1%  for  example  in  Case  C.  So  a  more  accurate 

analysis  would  be  to  assume  large  deflections  of  the  cables. 

The  boundary  conditions  assumed  for  the  cables,  how¬ 
ever,  are  the  most  representative  of  the  actual  situation. 
If  a  clamped  condition  were  to  be  assumed  at  the  input  link, 
the  system  would  be  stlffer  and  raise  K,  while  a  pinned  con¬ 
nection  at  the  bottom  would  lower  K,  yet  be  unrealistic  due 

to  the  firm  attachments.  Therefore,  in  comparison  with  the 


experimental  results  it  is  sbown  that  the  boundary  condi¬ 
tions  are  correct. 

From  the  comparison  of  the  experimental  and  analytical 
results  it  can  be  concluded  that  the  analytics  are  refined 
sufficiently  to  be  used  as  a  design  tool.  And,  as  discussed 
in  Chapter  7,  a  mechanism  can  be  designed  which  has  a  posi¬ 
tioning  error  of  less  than  one  arc  second. 


CHAPTER  6 


DAMPER 


6.1  Objectives 

As  previously  mentioned,  the  level  will  be  used  in  an 
environment  in  which  vibrations  are  present  due  to  the 
machine  tool  the  Instrument  is  monitoring  and  others  in  the 
surrounding  production  area.  These  vibrations  must  be  elim¬ 
inated  by  a  damper  because  the  leveling  mechanism  is  passive 
and  has  no  means  of  compensation.  The  characteristics 
looked  for  in  the  damper  were:  1)  the  ability  to  function 
equally  well  in  two  orthogonal  directions,  2)  remove  oscil¬ 
lations  above  1%  of  the  input  (100  arc  seconds  nominally) 
and  3)  not  interfere  with  the  function  of  the  level. 

These  constraints  lead  to  tne  cylindrical  damper  shown 
in  Chapter  1  Figure  1.4,  and  again  in  Figure  6.1.  Air  will 
be  used  as  the  viscous  medium  not  because  it  is  the  optimum, 
but  because  seals  are  not  required  and  hence  no  friction 
will  be  introduced  into  the  system  to  degrade  the  mechanisms 
positioning  accuracy.  The  analytics  developed  herein  will 
be  used  to  determine  the  feasibility  and  then  the  design  of 


•tv*e1 


a  cylindrical  air  damper  for  this  application 


6.2  Cylindrical  Pendulum  Damper  Macttanica 

Equation  6.1  [5],  the  Reynolds  Equation,  defines  the 
mechanics  of  a  fluid  between  walls  inclined  to  each  other 
and  moving  relative  to  each  other.  It  is  assumed  that  the 
space  between  the  walls  is  small  enough  to  suggest  that  the 
flow  is  laminar  and  fluid  velocity  perpendicular  to  the 
walls  is  negligible  compared  to  that  parallel  to  the  walls. 
These  conditions  prevail  approximately  if  the  ratio  between 
the  developed  flow  path  length  and  the  gap  is  better  than 
1000:1.  The  required  ratio  is  even  larger  if  the  fluid 
viscosity  is  very  low  as  is  the  case  with  air.  As  a  begin¬ 
ning,  approximations  based  on  laminar  flow  conditions  will 
be  assumed  even  though  the  length-gap  ratio  may  not  be  quite 
optimum. 


Reynolds  Equation: 


ac 


3 

h  dp 
6  tidy 


) 


dy 


f(x,y) 


f(x,y)  -  -  u|£  -  W~  +  2V  (6.1) 

See  list  of  symbols  for  Chapter  6. 

The  Reynolds  Equation  can  be  written  in  polar  coordi¬ 
nates  consistent  with  the  damper  geometry.  Furtner  assuming 
constant  viscosity  of  the  air 
IT]  Cameron,  A . ("Principles  of  Lubrication", 


J.  Wiley  Publ 


at  any  level  y  the  gap  h  can  be  written  as 


h  -  c  +  e(y)cos©  or  h  -  c  -e(y)cosP 
To  get  e  general  estlnate  of  the  mechanism  behavior.  It  Is 
advisable  to  non-d lmens lone  11 ze  Important  parameters.  For 
the  pendulum  damper,  the  following  transformations  are  con¬ 
venient. 


H(y) 


h(y) 

h 

o 


e(y) 


e(y) 

c 


^  (y-A)  -  h  (y-0) 
h  (y-0) 


0 


region 


h  (y-0)  -  fa  (y-A) 
h  (y-A) 


region 


Y 


L 

A 


In  a  pendulum  damper,  the  U  and  U  velocities  are  Inef¬ 
fectual  compared  to  the  squeeze  component  V  so  that  Equation 
6.2  can  be  simplified  and  with  the  non-dimensional  transfor¬ 
mation  becomes 


Mh3H)  / R »  2 

ae  va'  by 


-  -^~[-o  -o  +2 v j 


The  above  partial  differential  equation  can  not  be 


solved  in  a  closed  form  but  a  simplification  suggests  Itself 


T 


by  the  realities  of  tilting  dampers.  Because  b  varies  along 
y  It  Is  not  practical  to  use  long  damper  cylinders.  If,  on 

,  *  1  X 

the  other  hand,  short  cylinders  are  used  (  ~  <  ~  )  it  might 
be  assumed  that  pressure  gradients  In  the  circumferential 
direction  will  be  considerably  less  than  In  the  axial  dlrec- 
tlon.  If  that  is  the  case,  the  first  expression  on  the  left 
hand  side  of  the  Reynolds  Equation  might  be  ignored  reducing 
the  partial  differential  equation  to  the  ordinary  differen¬ 
tial  equation  6.3. 


where 
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the  squeeze  velocity  around  the  cylinder  is 
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A  similar  treatment  to  the  9  region  yields  an  analogous 
result  with  suction,  rather  than  positive  pressure,  being 
generated,  hence,  the  two  regions  reinforce  each  other  so 


that  the  total  fluid  force  becomes 


where 
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kd  Is  given  on  Figure  6.2  as  a  function  of  the  eccen- 

Jt 

tricity  ratio  e  and  length  ratios  —  .  At  e  ■  0  the  pendu- 

O  L  O 

lum  Is  exactly  centered  on  It's  cylindrical  housing.  At 
e  **  1.0  the  pendulum  Is  touching  the  wall.  The  minimum 
damping  force  occurs  at  e  ■  0 


There  were  several  approximations  made  to  expedite  the 
analysis.  While  the  results  are  not  exact,  the  errors  are 
probably  sufficient  for  design  purposes  and  they  do  show  the 
role  of  the  Important  parameters. 

6.3  Pendulum  with  Fluid  Damping 

As  a  convenient  but  reasonable  approximation  to  the 
actual  3  dimensional  dynamics,  a  2  dimensional  analysis  was 
developed.  A  sketch  of  the  damper  system  is  shown  in  Figure 
6.3. 


Taking  moments  about  the  pivot  0  (a  fixed  ins tantaneous 
center  of  curvature  rather  than  the  perturbated  small  move¬ 
ment  was  assumed): 


Figure  6.2  Damping  force  vs.  the  eccentricity 
ratio  e  as  a  function 
of  the  length  ratios  A/L. 
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or  for  small  motion, 
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The  undamped  natural  frequency  of  tne  pendulum  is: 
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The  damping  factor: 
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The  damping  force  is  described  as  F  - 


CDV  but 


Section  6.1  F.  is  also  described  as 
D 
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Substitution  of  this  into  Equation  6.9  yields 
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6 .4  Experiment 

Equation  6. 1.2  demonstrates  the  design  parameters  and 
their  role  for  acquiring  a  working  damping  coefficient. 
This  was  used  to  construct  the  experimental  apparatus  of 
Figure  6.4.  The  experiment  was  fairly  simple  in  that  the 
only  instrumentation  required  was  a  laser.  Tne  laser  was 
directed  to  a  mirror  attached  on  the  mechanism  and  then  on 
to  a  sheet  of  graph  paper  positioned  at  the  correct  height. 
By  virtue  of  the  distance  from  the  mirror  to  the  graph 
paper,  a  magnification  of  the  mirror  rotations  occur,  such 
that  up  the  laser  beam  was  displaced  about  6  cm.  for  a  one 
degree  swing  of  the  mirror.  Then  the  cusps  of  the  laser 
beam  were  marked  on  the  paper  as  the  time  was  recorded. 

The  experimental  damper  was  sized  as  shown  in  Figure 
6.5.  The  damping  ratio  was  then  determined  as  follows.  The 
clearance  (c)  is 

D-D 

c  “  - - -  ■  .0135 

Substituting  the  data  in  Figure  6.3  into  Equation  6.12 
yields 

2.56 

C  -  24  X  0.25  «  10-8  x  l.SJSt-^)3  «  2  *  «c 

or 

Q  -  0.11  Kd 


At  e  *  0,  K  ■  tt  (Figure6.2):  C  *  .014 


and  at  e  -  0.7;  X  -  -r—rr  -  .2;  K  -  0.85 
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It  X  is  increased  to  1"  and  every  other  property  Is 
unchanged  then 

„  .  13  ,2.84.2 

C  “  <.5)  *  ^.SO^  X  *lLKD 


at  e  •  0 


£  “  .14 


at  e  ■  0.7 


C  “  .93 


6.5  Experimental  Results 


The  first  experiment  was  conducted  with  X  ■  0.5".  The 
response,  Figure  6.4a,  indicates  £  “  0.17.  A  second  experi¬ 
ment  with  X  ■  1.0  yielded  the  response  of  Figure  6.4b  and 
£  ■  0.73.  The  experimental  set  up  did  not  allow  a  determi¬ 
nation  of  the  eccentricity  in  the  dormant  state,  further¬ 
more,  the  resolution  possible  was  coarser  than  £  -  0.3. 


Experimental  £  and  theoretical  £  are  0.17  and  0.094 
respectively,  from  the  first  test  (X  “  0.5").  Experimental 
and  theoretical  £  are  0.55  and  0.93,  respectively  from  the 
second  test  (X  ”  1.0").  Theoretical  values  were  based  on  a 
0.3  eccentricity  between  damper  sleeves  considering  that  the 
test  set  up  static  e  is  only  an  estimate  and  the  fact  that 
the  rapid  decay  of  oscillations  makes  it  difficult  to  meas¬ 
ure  the  second  and  third  decayed  amplitudes;  agreement  witn 
theory  seems  adequate  for  design  purposes. 


Equations  6.12  with  Figure  6.2  should  serve  as  adequate 


design  guides.  From  Equation  6.11  It  is  apparent  that  to 
maximize  damping  the  following  parameters  should  be  con- 
s idered . 

I 

1.  The  wetted  length  to  radial  clearance  ratio  —  is  very 
critical  coming  in  as  a  cubic  influence. 

2 

Lq 

2.  The  ratio  -  calls  for  setting  L  as  long  as  possible 

L1*5 

and  L  as  short  as  possible.  By  the  kinematics  of  the 
suspension  it  is  suggested  that  maximizing  L  may  serve 
both. 

3.  Both  the  mass  and  to  a  lesser  degree  the  radius  of 
gyration  should  also  be  minimized. 


Finally,  one  might  ask  how 
able.  In  other  words,  how  long 
subside  to  a  motion  1%  of  the 
decay  rate  is  approximated  by 

1 

n  “  —  In  100  ■  ' 


small  a  value  of  £  is  toler* 
will  it  take  the  pendulum  to 
initial  disturbance.  This 

In  100  .73 

"  C 


n  “  number  of  oscillations. 


If  the  period  of  the  pendulum  is  .6  secs  and  C  ■  . 1  the 
time  interval  is  4.4  seconds.  If  a  1  second  interval  is 
desirable  the  required  damping  factor  must  be  .44.  Hence  it 
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is  possible 
tling  time  of 
required  for 
damper  is  not 
both  nearly  3 


1 


to  design  an  air  damper  which 
1  second.  Also,  the  size 
this  example  is  reasonable, 
prohibitive  given  the  size  of 
73  inches. 
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CHAPTER  7 


CONCLUSIONS 

A  proposed  mechanism  to  level  optics  about  two  horizon¬ 
tal  axes  to  a  high  degree  of  accuracy  has  been  studied  to 
determine  the  structural  properties  that  influence  the 
operation.  The  mechanism  is  composed  of  two  ternary  links, 
one  being  the  input  with  the  other,  the  output,  being 
suspended  by  miniature  cables.  The  mechanism  was  designed 
to  simulate  a  spherical  pendulum  so  that  output  rotations 
would  be  less  than  1  arc  second  for  input  rotations  up  to 
100  arc  seconds. 

Initially,  the  3  dimensional  mechanism  was  generalized 
to  the  two  dimensional  case.  A  kinematic  analysis  was  per¬ 
formed  the  two  dimensional  model  to  highlight  the  effects  of 
link  lengths  and  tolerances  on  the  positioning  errors.  From 
this  analysis  it  was  learned  that  the  ratio  of  the  output 
link  length  to  the  input  link  length  was  an  important  param¬ 
eter.  This  ratio  must  be  maximized,  up  to  1000:1,  to  get  a 
desirable  performance  from  the  mechanism.  The  vertical 
position  of  the  center  of  mass  was  also  an  important  system 
parameter  and  should  lie  below  the  bottom  link  to  ensure 


minimum  kinematic  error.  A  surprising  result  from  this 
analysis  was  that  the  center  of  mass  need  not  be  centered  on 
the  output  link  due  to  the  assembly  procedure  of  the  mechan¬ 
ism. 

Upon  comparison  of  the  results  from  this  analysis  witn 
results  gathered  in  experiments,  it  was  found  that  the 
planar  model  predicted  lower  output  angles  for  two  of  the 
three  experimental  mechanisms.  It  was  hypothesized  that  the 
difference  was  due  to  rotations  of  the  output  link  which 
were  unaccounted  for  in  the  planar  model. 

A  spatial  kinematic  analysis  was  then  performed  which 
included  all  of  the  geometry  found  in  the  experimental 
mechanism.  While  the  analysis  was  far  longer  than  the 
planar,  a  small  amount  of  additional  information  was  gath¬ 
ered.  The  Important  result  was  that  the  actual  3  dimen¬ 
sional  mechanism  has  less  positioning  errors  than  Indicated 
by  the  planar  analysis  due  to  rotations  of  the  output  link 
about  the  vertical  axis  having  a  corrective  action  on  the 
rotations  about  the  horizontal  axes.  The  spatial  analysis 
checked  well  enough  with  the  planar  to  indicate  the  former 
could  be  done  with  little  loss  in  accuracy.  Since  the  spa¬ 
tial  analysis  indicates  lower  output  angles  than  the  previ¬ 
ous  analysis,  there  is  a  greater  discrepancy  with  the  exper¬ 
imental  data. 

The  effects  of  stiffness  in  the  mechanism  due  to  the 


bending  effects  of  the  support  cables  were  then  Investigated 


to  understand  the  difference  observed  between  the  analytical 
and  experimental  results.  The  support  cables  were  thought 
to  be  slender  beams  pinned  at  the  input  link  and  clamped  at 
the  output  link  on  the  two  dimensional  representation  of  the 
mechanism.  This  analytical  approach  proved  to  accurately 
predict  the  output  observed  in  the  experimental  work.  It 
was  found  that  the  effects  of  the  moment  of  Inertia  and 
modulus  of  elasticity  of  the  cables  were  stronger  parameters 
to  the  system,  for  certain  geometries,  than  the  kinematic 
properties.  It  was  shown  that  these  variables  increased  the 
output  angle  exponentially  such  that  the  minimum  modulus  and 
inertia  must  be  used  to  optimize  the  performance  of  the 
mechanism. 

An  air  damper  was  also  designed  to  minimize  distur¬ 
bances  of  the  optics  during  operation.  Two  cylinders,  one 
fixed  to  the  output  link  and  one  to  ground,  formed  a  two 
dimensional  damper  which  was  free  of  friction  while  dissi¬ 
pate  the  input  oscillations  effectively.  Again,  in  order  to 
establish  the  correct  geometry,  an  analysis  based  on  the 
Reynolds  Equation  assuming  constant  air  viscosity  was  per¬ 
formed.  Using  experimental  results  from  two  damper 
geometries,  it  was  shown  that  the  analysis  was  adequate  to 
characterize  the  system,  and  that  the  cylinder  length  should 
be  maximized  while  the  distance  from  the  input  link  to  the 
center  of  gravity  of  the  optics  should  also  be  maximized. 


Section  6.3  indicates  that  the  center  of  mass  should  he 


above  the  center  of  the  damping  force,  or  the  center  of 
the  wetted  Length  of  the  cylinder,  and  the  further  away 
the  better  for  maximized  damping.  From  Chapters  2  and 
4  one  would  like  the  center  of  mass  below  the  connec¬ 
tion  point  of  the  cables  to  the  output  link  to  minimize 
the  positioning  errors.  Hence,  the  center  of  mass 
should  lie  between  the  cable  connections  on  the  output 
link  and  the  center  of  the  wetted  length  to  satisfy 
both  the  settling  time  and  the  positioning  errors. 

The  centering  of  the  mass  on  the  ouput  link  was  shown 
to  be  unimportant  to  the  functioning  of  the  mechanism 
in  Chapters  2  and  4.  However,  the  miscenterlng  does 
effect  the  assembly  of  the  mechanism  to  the  damper. 
The  horizontal  shift  of  the  output  link  is  equal  to  the 
miscenterlng,  so  that  the  assembly  tolerance  of  the 
mechanism  to  the  damper  equals  the  tolerance  on  the 


Location  of  the  mass  center. 


APPENDIX  A 


Planar  Kinema  tic  Program  Lia  ting 

The  following  program  was  developed  to  determine  the 
equilibrium  position  of  a  four  bar  mechanism  with  a  mass 
attached  to  one  link.  The  program  calculates  the  wire 
lengths  (r2  &  r4)  required  to  make  the  output  link  (r3) 

level  Initially.  This  requires  the  user  to  supply  the 
lengths  of  the  input  (rl)  and  the  output  links,  an  initial 
rotation  (wlnt),  a  misscentering  of  the  optics  on  link  3  (u) 
and  the  height  of  the  mechanism  (h),  all  of  which  determine 
the  lengths  of  the  wires .  The  Instant  center  1/3  is  then 
positioned  directly  above  the  center  of  mass  on  the  output 
link.  A  non-linear  equation  solver  is  then  called  to  solve 
the  eight  constraint  equations  describing  the  location  of 
the  pin  joints,  points  A  and  B.  The  equation  solver  is  a 
routine  in  the  International  Mathematics  &  Statlsics  Library 
and  it  calls  an  external  subroutine  fen.  The  program  was 
written  in  Fortran  F77. 

program  main(input, output, tape5-lnput, tape6-output) 
c 

c  Planar  model  using  Newtonian  static  equilibrium, 
c  Length  of  wires,  links  2  and  4  are  calculated 
c  based  on  an  Initial  rotation  of  link  1  and  a 
c  miscenterlng  of  the  eg  on  link  3. 
c  Up-dated  on  Jan  2  1984 


dimension  par(l),wk(600),xo(8),xa(8) 

common/ da  ta/ rl ,  r2 , r3 ,r4 ,  r5 ,r6 , cx , cy , dx , dy 

real  gx,gy 

external  fen, fun 

n«8 

m- 150 

pi-acos ( - 1 .0 ) 
deg- 180 . /pi 

Variables  describing  link  lengths  are  read 
read(5,*)  r3 , r 1 , d 1 s t , ph , wint , u 
win-wlnt/deg 

Translate  link  1  until  the  IC  is  above  the  CG 
ax«-r3/2 .+u 
ay-0 

bx-r3/ 2 ,+u 
by-0 

cpx-ax*rl*(ph*cos(win)  +  bx*sin(win))/(bx*( ph+r 1*  s in( win) )-ph*ax) 
cpy-ph 

dpx-cpx+r l*cos (win) 
dpy-epy+r l*s ln( win) 
gx-0 .0 
gy-ay- d 1st 

Determine  the  lengths  of  links  2,  4,  5  &  6 
r2-sqrt((ax-cpx)**2.+(ay-cpy)**2.) 
r4-sqrt((bx-dpx)**2.+(by-dpy)**2.) 
r5-sqrt((gx-ax)**2.+(gy-ay)**2.) 
r6-sqrt(<gx-bx)**2.+(gy-by)**2.) 

Initial  estimate  of  independent  variables 
xo(  l)-ax 
xo ( 2 ) -ay 
xo(3)-bx 
xo ( 4 ) -by 
xo(5)-gx 
xo (6 )-gy 
xo(7)-l/2  . 
xo(8)-l/2. 

Output  indicating  initial  configuration 
wrlte(6,2) 

2  f orrna  t( lhl ) 
write (6, 10)  rl 

10  f orma t( 6x , 13hor if ace  dia  -,f7.4) 
write(6,20)  r2,r4 

20  f orma t( 6x  ,  1 7hleng th  of  wires  -,f7.4) 
write (6, 30)  r3 

30  f orma t( 6x , 8hs tance  -,t7.4) 
write (6, 35)  dist 

35  f orma t( 6x , 13 hp la te  to  cg.-,f7.4) 
write (6, 36)  ph 

36  f orma t( 6x ,8hheigh t  -,f7.4) 
write(6,38)  wint 

38  f orma t ( 6x  ,  2 Ihinl tla 1  angular  error, f7. 4) 


wr 1 te ( 6 , 39  )  u 

39  £ormat( 6x ,9hcg .  error, £7. 4) 
wr  i  te  (  6 ,40  ) 

40  format(/// , 8x , 7  h  the  ta  g,6x,7htheta  m ,6x , 7 h the ta  2 
6x,7htheta  3,5x,9hd  theta  m) 

read  thd  number  o£  Input  states  to  be 
solved,  1.  Read  the  Input  state  tlnt(deg) 
read ( 5,3)  1 
5  £ormat(l2) 
do  100  j-1,1 
read(5,*)  tint 
tin*tint/deg 
cx»cpx+rl*cos ( t in ) / 2 . 
cy«cpy-rl*sin(tin)/2. 
dx»dpx-rl*cos(tin)/2. 
dy«dpy+rl*sin(tin)/2. 

Non-linear  equation  solver 
call  2scnt(fcn,10,n,m,par,xo,fnorm,wk,ier) 

Calculate  the  rotation  of  link.  3 
from  points  ay  and  by. 
ay**xo  (2) 
by«*xo  (  4  ) 

tho*asin((by-ay)/r3) 

tout*tho*deg*3600. 

Calculate  the  input  output  ratio  rat 
1  £ ( tint . ne . 0 . 0 )  then 
rat»tho*deg/tint 
else 
ra  t«*0 
endi£ 

write(6,60)  tint,tout,rat,xo(7) ,xo(8) 
b 0  forma t ( / , 3x , 5  £  13 . 5 ) 

100  continue 

s  top 
end 

Subroutine  containing  constraint  equations 
subroutine  fen  (x,£,n,par) 

common /da ta/ ri ,r2,r3,r4,r5,r6,cx,cy,dx,dy 

dimension  x(n) , f (n) , par( 1) 

real  gx,gy 

ax-x ( 1 ) 

ay-x<2) 

bx«x(3) 

by-x(4) 

gx-x<5) 

gy-x(6) 

f  ca**x  (  7 ) 

£  d  b-x ( 8) 


f  (l)-(ax-bx)**2.  +  ( ay- by ) **2 .  -  r3**2. 
f ( 2 )- (ax-cx) **2 .  +  (ay-cy)**2.  -  r2**2. 
f ( 3 )■ ( bx-dx) **2 .  +  (by-dy)**2.  -  r4**2. 
f (4)- (ax-gx)**2 .  +  (ay-gy)**2.  -  r5**2. 
f (5)-(bx-gx)**2.  +  (by-gy)**2.  -  r6**2. 
f  (6)*( cx-ax)*£ca/r2  +  (dx-bx)*f db/r4 
f ( 7 ) » ( cy-ay) *f ca/r2  +  ( dy- by ) * f d b/ r4  -1. 
f ( 8 )■ ( ( bx-ax) * ( dy- by )  -  (  by-ay )*( dx- bx ))* f d b/ r4 
-(gx-ax) 
re  turn 
end 
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APPENDIX  B 


Spatial  Kinematic  Program 

Tula  program  was  developed  to  perform  a  spatial 
kinematic  position  solution  for  a  five  bar  mechanism  with 
spherical  joints  and  a  mass  fixed  to  one  link.  The  Input  to 
the  program  are  the  lengths  of  the  input  links  (rl,  rZ ,  r3), 
the  height  of  the  mechanism  (ph),  the  lengths  of  the  output 
links  (r4,  r5 ,  r6)t  the  miscentering  of  the  eg  along  the  X 
axis  (u)  and  the  vertical  distance  from  the  output  plate  to 
the  eg  (dlst).  The  program  calculates  the  wire  lengths 
based  on  the  above  variables  and  an  Initial  rotation  of  the 
top  plate. 

A  subroutine  (rotate)  calculates  the  final  position  of 
the  top  spherical  joints  via  a  rotation  matrix.  The  output 
rotations  are  calculated  based  on  Rodrigue's  formula  and  the 
normal  to  points  d,  e,  and  f.  As  before,  a  non-linear  equa¬ 
tion  solver  is  used  to  find  the  equilibrium  of  the  mechanism 
and  requires  subroutine  fen, 

program  main(lnput,output,tape5"lnput,tape6aoutput) 
real  nx,ny,ns,h,hp,knx, kny , knz , kn 

dimension  xo(15),wk(1104),par(l),xp(15),h(3),hp(3) 


o  o 


!  ,ztp(3),ytp(3),xtp(3),xt(3),yt(3),*t(3) 

eoaaoo/dttt/rl ,r2,r3,r4,r5,r6,r7,r8,r9,rl0,rll,rl2 
!  ,xt,yt,xt 
axtaraal  fen 
e 
e 

e  xd-x(l)  xe«x(4)  xf»x(7)  xeg-x(lO)  fda-x(l3) 

e  yd-x(2)  ye»x(5)  yf-x(8)  ycg-x(il)  feb-x(14) 

c  zd-x(3)  *e«x(6)  zf-x(9)  zcg-x(l2)  ffc«x(15) 


n-15 

•■150 

c 

pl>icos ( -1 .0 ) 
deg-180 ./pi 
e 

c  initial  conditions,  calculate  the  geometric  center 

c  of  the  top  plate  and  then  points  a,  b,  and  c. 

c 

read ( 5,*)  rl,r2,r3 
read(5,*)  r4,r5,r6 
read(5,*)  dist,ph,u 

Phil2-acos((rl**2.+r2**2.-r3**2.)/(2*rl*r2)) 
Phil3-acos<(rl**2.+r3**2.-r2**2.)/(2*rl*r3) ) 
phi23-pl-phll2-phil3 

phl45-acos((r4**2.+r5**2.-r6**2.)/(2*r4*r5) ) 
phi46-acos( (r4**2 .+r6**2 .-r5**2 . ) / ( 2*r4*r6 ) ) 
phi56»pi-phl45-phl46 
c 

tp>pi/2-phl23 
rp»rl/(2*cos(tp)) 
t b«pl/2-ph!56 
rb-r4/(2*cos(tb)) 
c 

xtp(l)--rp*cos(tp) 
y tp( 1)-0.0 
ztp(l)*rp*sln(tp) 
xtp(2)-0.0 
ytp(2)-0.0 
z  tp( 2 ) *-rp 
xtp(3)-rp*cos(tp) 
y tp(3)-0 .0 
ztp(3)*rp*sln(tp) 

c  Initial  estimate  for  points  d,  e,  f  and  g,  also  forces 

xo(l)»-rb*cos(tb) 
xo(2)»ph 
xo(3)-rb*sln( tb) 
xo (4)-Q  .0 
xo(5)-ph 
xo(6)»-rb 
xo(7 )-rb*cos( tb) 
xo(8)«ph 


xo(9)-rb*sln(tb) 
xo(10)-u 
xo( ll)-ph+dlst 
xo(12)-0 .0 
xo(13)-0. 33333 
xo(14)-0. 33333 
xo(i5)-0. 33333 

Normal  to  lower  piece  Initially 
h(l)-xo(l) 
n( 2)-xo(2)+i . 
h( 3)-xo( 3) 

Score  for  initial  position 
do  11,1*1,15 
11  xp(i)«xo(l) 

Initial  rotation  ot  top  plate 
read(5,*)  wxlnt,wzint 

call  roCate<xtp,ytp,ztp,xt,yt,zt,wxint,wzlnt) 
do  90  1-1,3 
xtp(l)— xt(i) 
ytp(i)-yt(l) 
ztp( 1 )  —  z t( 1) 

90  continue 

Calculate  wire  lengths  and  CG  constraint  links 

r7-sqrt((xod)-xtpd))**2.+(xo(2)-ytpd))**2.  +  (xo(3)-ztpd))**2.) 
r8-sqrt((xo(4)-xtp(2))**2.+(xo(5)-ytp(2))**2.  +  (xo(6)-ztpd))**2.) 
r9-sqrt( (xo(7 )-xtp(3) )**2.+(xo(8)-y tp(3) )**2.+(xo(9)-ztp(3))**2.) 
rl0-sqrt((xo(10)-xo(l))**2.+(xo(ll)-xo(2))**2.+(xo(12)-xo(3)) 

!  **2 . ) 

iU«sqrt((xo{10)-xo(4))**2.+(xo(U)-xo(5))**2.+(xo(12)>xo(6)) 

!  **2 .  ) 

rl2-sqrt((xo(10)-xo(7))**2.+(xo(ll)-xo(8))**2.+(xo(l2)-xo(9)) 

!  **2  . ) 

Input  rotation  of  top  plate,  do  (1)  times 
write(6,170) 
wr lte ( 6 , 17  5) 
r ead ( 5,10)  1 
10  format(l2) 
do  350  j-1,1 
read(5,*)  omegx,omegz 

call  rotate(xtp,ytp,ztp,xt,yt, z  t , omegx , omegz) 

Call  equation  solver 

call  zscnt  ( f cn ,9 , n ,m , par , xo , t norm , wk , ler ) 

Rodrigue's  formula 

xha-(xo(8)-xo(2) ) * ( xo ( 6 ) -xo ( 3 ) )-(xo(5)-xo(2))*(xo(9)-xo(3)> 

y  ha-  (  xo ( 4 ) >xo ( l))*(xo(9)>xo(3))-(xo(7)-xo(l))*(xo(6)-xo(3)> 

zha-(xo(7 ) -xo  ( 1) )*(xo(5)-xo(2) )-(xo(4)-xo( 1) )*(xo(8)-xo(2) ) 

tl-sqr t(xha**2 .+yha**2 .+zha**2 . ) 

hp(l)-xo(l)+xha/ tl 

hp( 2 )-xo( 2 )+yna/ tl 

hp(3)-xo(3)+zha/tl 

knx-(xo(2)-xp(2))*(hp(3)-h(3))-(hp(2)-h(2))*(xo(3)-xp(3)) 

kny-(hpd)-hd))*(xo(3)-xp(3))-(xod)-xp(l))*(hp(3)-h(3)) 


knz-(xo( 1) -xp( 1) )*(bp(2)-b(2) )-(bp( i)-b(i) )*(xo(2)-xp(2) ) 

kn-sqr t(knx**2 .+kay**2 .+knz**2.) 

nx-knx/kn 

ayakny/kn 

az-kaa/kn 

sl-(b(l)-xp(l))*(bp(l)~xo(l))+(b(2)-xp(2))*(bp(2)>xo(2))+(b(3) 

xp(3))*(bp(3)-xo(3)) 

s2«(b(l)-xp(l) )*nx+( h( 2 ) -xp( 2 ) ) *ay+( b( 3) -xp( 3) )  *nz 

s3«(b(2)-xp(2))*az-(b(3)-xp(3))*ay 

s4-(b(3)-xp(3) )*ax-(ta(l)-xp( 1) )*nz 

s5«(h(l)-xp( l))*ny-(ta(2)-xp(2) )*ax 

a6«s3**2.+s4**2.+a5**2. 

b«ta>aco(( (sl-(s2**2.))/s6) 

thetax-beta*ox*deg*3600. 

thetay«beta*ny*deg*3600. 

thetaz-beta*nz*deg*3600. 

thdm-sqr t(  thetax**2.+thetaz**2.) 

tbln-sqrt(omegx**2.+omegz**2.) 

Calculate  laput  output  ratio  rat 
if( thia.ae .0 .0 )  then 
rat-thdm/(  tb in* 3 6 00  .) 
else 

rat-0.0 

Create  new  aormal  if  tbe  first  one  was  superseded 
h ( l ) -xo ( 1 ) 
h( 2 )-xo( 2 )+l . 
b(3)-xo(3) 
do  305,  1-1,15 
•  05  xp(l)-xo(l) 
eadlf 

write (6,174)  thdm , ooegx , omegz , ra  t , f nor a , tfte  tay 
50  continue 

ou  tpu  t 

write (6, 130)  rl,r2,r3 
write (6 , 160)  r4,r5,r6,pb 
write (6, 190)  r7,r8,r9 
write(6 , 135) 

write (6, 140)  (x tp ( 1 ) , 1- 1 , 3 ) 
wr 1 te ( 6 , 140 )  (y tp( i ) , 1- 1 , 3 ) 
wrlte(6,140)  ( z tp( 1 ) , 1- 1 , 3 ) 
write(6 ,145) 

wr 1 te( 6 , 140)  (x t( i ) , i- l , 3 ) 
write (6, 140)  (y t ( i ) , i- 1 , 3 ) 
write (6, 140)  ( z t( 1 ) , 1- 1 , 3 ) 
write(6 ,165) 

write (6, 140)  xo(l),xo<4),xo(7),xo(10) 
write(6,140)  xo(2),xo(5),xo(8),xo(ll) 
write (6, 140)  xo(3),xo(6),xo(9),xo(12) 
wrlte(6,180)  wxlnt,wzlnt 

.30  f oraa t( lbl5x , 'upper  plate  diBeas ions : ' , / iOx ,' r 1«" , f 10 .4 , 

1  3x,'r2-',fl0.4,3x,'r3-',f 10.4) 


.45  forma t(/  ,  22x ,'a"  ,  14x , "  b" ,  14x,"c") 

.40  format( L0x,4f 15 .5) 

.35  format(/,21x,"ap",13x,"bp",l3x,"cp") 

.50  format(/ ,5x , 'input  rotation  was:  (deg . ) " , / lOx , 2fl0 . 2) 
.80  format(/,5x, "lower  plate  dimens  ions :",/ lOx  ,"r4-" 

!  ,f7.3,3x,"r5-",f7 .  3 , 3x  , " r6-" , f 7 .3 , 3x,"h»" ,f7 .3) 

L  65  forma  t( /2 lx ,"d",14x,'e",14x,"f",14x,"cg") 

L  7  0  f  orma  t  ( 16x  ,  "  tiie  ta  out", 

1  ,6x  ," omega x" ,6x ," omegas" , 6x ," ra t" ,6x ,"f norm" ) 

L75  format(18x,"  see.  ",6x,"  sec.  ",6x,"  see.  ") 

L  7  4  f orma t ( / , lOx ,6f 12 . 5 ) 

L80  format(/ ,5x , "initial  angular  of t se t" , / , lOx , 2£ 10 . 2 ) 

L  90  format(/,5x,'wire  lengths" , / , lOx ,"r7«" , f 7 . 3 ,3x , 

I  "r8-",f7.3,3x,"r9-',f7.3) 
s  top 
end 

Subroutine  used  to  rotate  top  plate 

subroutine  rotate(xl,yi,zi,xf,yf,zf,omegax, omega z ) 
d imension  xl(3),yl(3),zl(3),xf(3),y£(3),zf(3) 
deg«180./(acos(-1.0)) 

1 f ( ( omega x+omegaz ) .It. 0.00001)  then 
do  90  1-1,3 
xf (i)-xi(i) 
yf (i)-yi(i) 
zf ( i)-zi(i) 

90  continue 
else 

omegx-omegax/deg 
omegz-omegaz/deg 
omega-sqrt(omegx**2 .+omegz**2.) 
ux-omegx/ omega 
uz-omegz/ omega 
do  10  j-1,3 

xf(J)-xl(j)*((ux**2.)*(l-cos(omega))4-cos(omega)> 

!  +zi( j )*(ux*uz*( l-cos( omega )))-yi(j )*s in( omega ) *uz 

yf(j)-xl(j)*uz*sln(omega)-zi(j )* ( ux*s in ( omega ) ) 

1  +yi(j )*cos( omega) 

zf(j)-xl(j)*ux*uz*(l-cos(omega))+zi(j)*((uz**2.) 

!  *(l-cos( omega ))+cos(omega))+yl(j )*ux*s in ( omega) 

10  continue 
endlf 
re  turn 
end 

Subroutine  used  by  zscnt 

subroutine  fen  (x,f,n,par) 
dimension  x( n) , f ( n ) , par ( 1 ) 

!  ,  x  t  (  3 )  ,  y  t (  3  )  ,  z  t(  3 ) 

common/ data  /  rl  ,r2,r3,r4,r5,r6,r7,r8,r9,r!0,rll,rl2 
!  ,xt,y t, zt 


f (l)-((xt(i)-x(l))*x<13)/r7)+((xt(2>-x(4))*x(14)/r8)+((xt(3) 

!  -x(7) )*x( 15)/r9) 

f (2)-((yt(i)-x(2))*x(13)/r7)+((y t( 2 ) -x( 5 ) )*x( 14) /r8 )+( (y t( 3) 

!  -x(8))*x(15)/r9)-1.0 

f(3)-((zt(l)-x(3))*x(13)/r7)+((zt(2)-x(6))*x(14)/r8)+((zt(3) 

!  -x(9))*x(!5)/r9) 

f(4)«(((x(2)-x(ll))*(zt(l)-x(3))-(x(3)-x(i2))*(yt(l)-x(2))) 

!  *x(13)/r7) 

!  +(((x(5)-x(ll))*(zt(2)-x(6))-(x(6)-x(i2))*(yt(2)-x(5))) 

!  *x(14)/r8) 

!  +(<(x(8)-x(ll))*(zt(3)-x(9))-(x(9)-x(12))*(yt(3)-x(8))) 

!  *x(15)/r9) 

f (5)-(((x(3)-x(12))*(xt(l)-x(l))-(x(l)-x(lO))*(zt(i)-x(3))) 

!  *x(13)/r7) 

i  +(((x(6)-x(I2))*(xt(2)-x(4))-(x(4)-x(iO))*(zt(2)-x(6))) 

!  *x(14)/r8) 

!  +  <((x(9)-x(12))*<xt(3)-x(7))-(x(7)-x(10))*Ut(3)-x(9))) 

!  *x( 15) / r9 ) 

£(6)-(((x(l)-x(10>)*(yt(t)-x(2))-(x(2)-x(ll))*(xt(l)-x(l))) 

!  *x( 13)/r7) 

i  +(((x(4)-x(10))*(yt(2)-x(5))-(x(5)-x(ll))*(xt(2)-x(4))) 

!  *x( 14 ) / r8 ) 

!  +(((x(7)-x(10))*(yt(3)-x(8))-(x(8)-x(ll))*(xt(3)-x(7))) 

!  *x(l5)/r9) 

t <7)-<x<l)-xt<l))**2.+(x(2)-yt<l))**2.+(x(3)-zt(i))**2.-(r7**2.) 
f (8)»(x(4)-xt(2))**2.+(x(5)-yt(2) ) ** 2 .+(x( 6 ) - z t (2) )**2 . - ( r8**2 . ) 
f (9)»(x(7)-xt(3))**2.+(x(8)-yt(3))**2.+(x(9)-zt(3))**2.-(r9**2.) 

f (10)-(x(4)-x(l))**2.+(x(5)-x(2))**2.+(x(6)-x(3))**2.-(r5**2.) 
f  (U)-(x(7)-x(l))**2.  +  <x(8)-x(2))**2.+  (x(9)-x(3))**2.-(r4**2.) 
f (12)-(x(7)-x(4))**2.+(x(8)-x(5))**2.+(x(9)-x(6))**2.-(r6**2.) 

£<13)«<x(10)-x(l))**2.  +  <x(U)-x(2))**2.+  (x(12)-x(3))**2.-ri0**2. 
t(14)-(x(10)-x(4))**2.+(x(ll)-xC5))**2.+(x(12)-x(6))**2.-rll**2. 
f  (15)-(x(10)-x(7))**2.  +  (x(H)-x(8))**2.  +  (x(12)-x(9))**2.-rl2**2. 


APPENDIX  c 


Planar  Stlf f nan  Program 


This  program  calculates  the  equilibrium  position  of  a 
four  bar  mechanism  with  flexible  members  on  links  2  &  4  and 
a  mass  attached  to  link  3.  The  user  reads  In  the  values  for 
mg,  El,  the  lengths  of  links  1,  2,  3,  and  4  (rl,  r2 ,  r3 ,  & 
r4)  as  well  as  the  Input  angle  (thtl  deg.)  for  each  Input 
state.  As  with  the  earlier  programs  a  subroutine  Is  called 
to  solve  the  9  constraint  equations  for  the  position  solu¬ 
tion.  The  program  is  written  In  Fortran  F77. 


c 

c 

c 

c 

c 

c 

e 

c 

c 


Planar  analysis  with  the  wires  modeled 
as  small  deflected  beams.  3/16/84 
The  boundary  conditions  are  pinned  at  the 
top  and  clamped  at  the  bottom. 


xl*tho  x4-va  x7*mb 
x2-ha  x5-vb  xS-da 
x3**hb  x6“oa  x9*db 


program  maln(lnput,output,tape5-lnput,tape6»output) 
dimension  x( 9 ) , par ( 1 ) , wk( 1000 ) 
real  ll ,  la , lb , lo  ,mg 
integer  i,j 

common  u,d,el,thl,tha,thb,deg,ll,la,lb,lo,mg,j 
external  f cn , cosh , s lnh 


pl-acos(-l  .0) 
deg-180  ./ pi 
ltnax-125 
n-9 


c 


read ( 5  ,  * )  mg,el,d,u 


el-4.94e-4 

Initial  iitlnate  tor  aquation  aolvar 
x(l)-0.0 
x(2)  — 0.0847 
x(3)-0 .0847 
x( 4) -mg/ 2 . 
x(5)-mg/2  . 
x(6)-0.0 
x(7)-0.0 
x ( 8 )-0 .0 
x(9)-0.0 

User  inputs,  do  (1)  times 
read (5, 2)  1 
forma  t ( 12 ) 
do  40  1-1,1 

r ead ( 5  ,  *  )  ll  ,  la , lb , lo , thtl 
thl-thtl/deg 

Calculate  initial  angles 
lf(l.ne.l)  goto  5 
xa- (  lo- 1 i ) / 2 . 
tha-asin(xa/ la) 
ttib-as  ln(xa/  lb) 

5  continue 

Call  equation  solver 

call  zsent(fcn,13 ,n  ,  1 tmax , par , x , f norm , wk , ler) 
Calculate  output 
da-x(8) 
d  b-x( 9) 
pba- the- thi 
phb-thb+thl 
tot-x( 1 ) *deg 

Calculate  input  output  ratio 
rat-0.0 

if  ( thtl . le .0 .0 )  goto  11 
rat-tot/  tntl 
11  continue 

xb-li*cos( tbi)+lb*sin(phb)-do*cos(pnb) 
yb-li*ain(thi)-lb*cos(phb)-db*aln(pbb) 
xa--la*sln(pba)-da*cos(pba) 
ya--la*cos(pha)+da*sln(pha) 

Output  critical  variables 
wrlte(6,20)  th tl , to t*3600 . , ra t 
write (6, 20)  li , lo , la , 1 b , e 1 
10  format( 10x,4el3 .4) 
vrite( 6 ,30 )  (x(j),j-l,7) 
wrlte( 6,30)  (x(j),J-8,9) 

write (6, 30)  xa,ya,xb,yb,pba*deg,pbb*deg 
30  formet(5x,7f 12 .5) 

20  format(8x,7f 10.4) 
j-0 

40  continue 


Subroutine  used  by  zscnt 

subroutine  fen  (x,f,n,per) 
dimension  x( n) , f (n) , per ( 1 ) 

common  u  ,  d  ,  el ,  thl ,  tbs  ,thb,deg,li,la,ib,lo,mg,j 
reel  11 , le , lb , lo ,ma ,ob ,mg , be , bb 
external  cosb,sinb 
Integer  i,J 

tho«x( 1 ) 
be»x( 2) 
hb-x( 3) 
va«x(4) 
vb-x( 5) 
oe*x( 6) 
mb«x( 7 ) 
da«x(8) 
db*x(9) 

pha- ttia-  thl 
phb-thb+thl 
bl«sqrt( ve/ei) 
b2*sqrt( vb/el) 
cblle>cosb(bl*le) 
sblla-slnh(bl*la) 
cb21b"cosb(b2*lb) 
sb21b«slnh(b2*lb) 

f(l)«,ha*cos(pha)  +  va*sin(phe)+hb*cos(phb)-vb*sin(phb) 
f (2)«va*co* ( p be) -be *s in ( phe )+ho*s m (phb)+vb*cos ( p hb) 
f ( 3)«-me-mb+vb*lo*cos( phb- tho)+hb*io*sin(phb- tho) 
-mg*((lo-u)*cos(tho)/2.+d*sin(tho)) 

£(4)»ha*sblla/( bl* c bl la )+ma 
f(5)-hb*sb21b/(b2*cb2lb)+rab 

f(6)»-ha*sbila/(va*bl*cblia)+( tho+ha/ ve) *  la -da 
f(7)«-hb*sb21b/(vb*b2*cb2lb)+( tho+bb/ vb )* lb-d  b 

xa»-la*sln(ptia)-da*cos(pha) 

ya»-la*cos(pha)+da*8in(pha) 

xb*ll*cos(thi)+lb*sin(phb)-db*cos(phb) 

yb"li*sln(thi)-lb*cos(phb)-db*sin(phb) 

f(8)"(xb-xa)**2.+(yb-ya)**2.-lo*»2. 

£(9)"(yb-ye)/(xb-xa)-tan(tbo) 

return 

end 


Evaluate  the  slab  and  cosb  m  functions 


function  cosb(e) 


real 

e 

eoib" 

(exp(a)+exp(-a) )/2 

re  turn 

end 

f  unet 

ion  alnh(a) 

reel 

e 

a  inh« 

(exp(a)-exp(-e) )/2 

return 
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